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Symbols 


U noise 

measurement noise weighting 

a 

bandwidth of elcvon dynamics. 


matrix 


30 rad/sec 

Wpl • ] 

performance weighting matrix for 

b 

bandwidth of fuel flow rate dynam- 
ics, 100 rad /sec 

Wa 

given variable 

uncertainty weighting matrix 

C 

set of complex numbers 

w 

uncertainty matrix input variable 

D 

positive definite Hermit ian matrix 
exogenous inputs, (i.e., noise. 

Wj 

unity magnitude white Gaussian 

d 


noise with zero mean 

commands, turbulence) 

d' 

generalized exogenous input, 

^ noise 

state measurement, noise variable 


d'=[zd] 

X 

state longitudinal variable 4 

V 

error/ performance 

y 

sensed variable 4 

(■' 

generalized output performance 

' y =M 

z 

uncertainty matrix output variable 4 

Fi(P.K) 

closed-loop response between P 

a 

iiHTtial angle of attack. de 4 g 


and d f 

a 

peak magnitude of frequency 

F„(G, A) 

perturbed closed-loo]) response 4 


response 


between e and d 

A 

uncertainty matrix 

F u , Fa- 

longitudinal and vertical Drydem 
turbulence filters 

At 

time 4 increment 

f 

stoichiometric ratio for hydrogen, 

A 

general se 4 t of complex perturbations 

0.029 


not restricted by bounds 

G{. s) 

augmented system plant that 

Ab(\ A brii f 

uncertainty in control effee*t.ors 


contains /\(.s)and P(s) 

Banach space consisting of funct ions 

symmetric elcvon, elt 4 g 

be 

H 0 c 

be r , briij r 

e 4 ont TollcT-commanded control 

F(s) which are analytic in Re s > 0. 


take 4 values in C" x/ " , and are 
bounded in Re s > 0 

inputs 


h 

altitude 4 , ft 


effective control inputs 

In 

n x n identity matrix 

bc r , bjllj f. 

uncertainty matrix inputs in general 
structure 

K(.s) 

controller 



M 

Mach number 

V 

fuel equivalence ratio 

din 

air mass flow rate 4 , slugs/sec 

0 

pitch angle, ele 4 g 


fuel mass flow rate, slugs/sec* 

/' 

struct ureel singular value 

P(.H) 

generalized plant structure 

<7 

maximum singular value 4 

Q 

unitary matrix 

Abbreviations: 


<1 

pitch rate, de 4 g/se 4 e 4 

DOF 

degree of freenlom 

Re 

,S 

real part- 

variable e)f Laplace multiplier 

iff 

if and only if 

a 

control effectors variable 

LTI 

linear time 4 invariant 

V 

velocity, ft /sec 

SSTO 

single stage to e>rbit 


V 
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Notation: 

Hp 

[V<h] r 

[V,h] <v 

II ' Ik 


performance- weighted variables 
cominanded velocity and altitude 

performance-weighted velocity and 
altitude error 

infinity norm 


|| • ||^ maximum structured singular value 

over all frequency u> 

sup least upper bound over all frequency uj 

ijJ 

A dot over a symbol indicates differentiation with 
respect to time. 
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Introduction 

Single- stage- to- or bit (SSTO) air-breathing vehi- 
cles present significant challenges in many techno- 
logical areas and especially in flight control. These 
vehicles traverse a broader flight envelope than any 
aircraft flown to date and performance is critical 
throughout the entire flight regime to achieve the 
mission objective. Large variations in vehicle static 
and dynamic characteristics and in mass properties, 
such as significant movement of the aerodynamic 
center of pressure, result in continuously changing 
static stability margins throughout the flight enve- 
lope. (See ref. 1.) Furthermore, an additional source 
of uncertainty arises from the accuracy of mathe- 
matical dynamic models used to describe the vehicle 
in control system design. These challenges and the 
limited availability of empirical data above a Mach 
number of 8 in aerodynamics, propulsion, aeroelastic- 
ity, and heating as well as limited knowledge of their 
combined effects on vehicle mission performance dic- 
tate the need for a. robust yet performance-oriented 
control system. 

The airframe/propulsion interactions, possibly 
the most complex of any vehicle, an? critically im- 
portant to a successful hypersonic vehicle mission. 
The high sensitivity of the air-breathing propulsion 
system performance to changes in angle of attack and 
to dynamic pressure has been identified (ref. 2) and 
confirmed. (See ref. 1.) Furthermore, atmospheric 
turbulence and large density variations at a high al- 
titude and at a high Mach number (ref. 3) introduce 
other significant sources of uncertainty in the perfor- 
mance of an air-breathing propulsion system for the 
control system. 

In addition, such parameter uncertainties as 
propulsive efficiency, drag, and weight have major ef- 
fects on performance margins as the vehicle reaches 
orbital speed. (Sec ref. 4.) The significantly detri- 
mental effect of control surface deflection-induced 
drag on fuel consumption to orbit (ref. 5) provides 
another compelling reason for control system opti- 
mization in hypersonic vehicles. 

These issues and their impact on control system 
development have been recognized by numerous re- 
searchers. (Sec refs. 1, 6, and 7.) The control work 
in this area has primarily addressed an air-breathing 
hypersonic cruiser (ref. 7), which assumes equilib- 
rium steady-state flight with changes in the coeffi- 
cients of the equations of motion stemming from poor 
model description rather than from changing flight 
parameters due to acceleration. Furthermore, the 
control laws developed for the ascent phase have ei- 
ther disregarded the impact of angle-of-attack varia- 


tions on air-breathing propulsion performance (ref. 7) 
or, while addressing tracking and atmospheric turbu- 
lence issues, did not explicitly consider performance 
robustness. (See ref. 8.) 

Recent application of modern robust control the- 
ory to Space Shuttle (ref. 9) and fighter airplane 
(ref. 10) flight control problems demonstrated poten- 
tial benefits to the above challenges. The ob jective of 
this research is twofold: to assess applicability and to 
exploit the capability of modern multivariable robust 
control theory to explicitly deal with performance 
and with the uncertainty arising from changing flight 
conditions and vehicle characteristic s. The problem 
is formulated to deal with the challenge's associated 
with an SSTO vehicle and air-breathing propulsion 
system. A structured uncertainty model, represent- 
ing parametric variations as actuator uncertainty, is 
used to illustrate* the distinction between two modern 
design procedures, Hoc fttid // synthesis. 

In the first section of this paper is a brief re- 
view of robustness measures. The issue of possi- 
ble* conservative solutions to some practical prob- 
lems is discussed and the structuml singular value 
is introduced. The theoretical basis for the two syn- 
thesis procedures is also provided. In the next sec- 
tion is the derivation of a five-state longitudinal lin- 
ear model description from an ascent trajectory of a 
conical accelerator vehicle. After that is the prob- 
lem description of the uncertainty model, the Hoc 
weighting function selection, and the explicit inclu- 
sion of stochastic atmospheric t urbulence in a con- 
troller design. The last section deals with the appli- 
cation of Hoc and p controller synthesis and analysis 
techniques. The conclusions derived from this study 
follow. Appendixes A E provide further details of 
the hypersonic vehicle model and controller design. 
Included in the appendixes are an updated scramjct 
propusion model, the derivation of linear models from 
points along a continuous-acceleration flight, trajec- 
tory, an atmospheric model implementation in the 
design and simulation, and the numerical values of 
the original plant and the derived controllers. 

Theoretical Review of // 

This section provides essential theorems for ro- 
bustness and performance analysis in a system with 
uncertainty. Far more detailed and rigorous discus- 
sion is presented in references 11 13. 

Analysis methods based on singular values have 
been successful in providing multiloop extensions for 
classical single-loop techniques. (See ref. 11.) How- 
ever, these methods are limited to providing nec- 
essary and sufficient conditions for robust stability 



for systems with unstructured uncertainty, defined 
as norm-bounded but otherwise unknown perturba- 
tions. Consider, for example, the standard problem 
of analyzing a feedback system with simultaneous 
multiplicative uncertainty at the plant input and out- 
put signals. To apply singular value techniques, both 
perturbations must be reflected to a single location in 
the feedback loop, thus immediately inducing conser- 
vatism. However, because the combination of linear 
transformations is linear, any uncertainty occurring 
at several different locations in the feedback loop can 
be rearranged as a single block-diagonal perturbation 
in a larger feedback loop. In other words, even un- 
structured uncertainty at the loop component level 
becomes highly structured at the system level. (From 
MUS\ N: Robust Multivariable Control Short Course 
and Software, available from John C. Doyle, Andrew 
Packard, and Gary Balas, MUSYN, Inc.. Minneapo- 
lis, Minnesota.) 

The general framework for the robust control 
problem is introduced in figure 1(a). Any linear 
combination of inputs, outputs, commands, pertur- 
bations, and controller can be arranged into the form 
in the diagram. (See ref. 12.) Furthermore, the ex- 
ogenous inputs d, the perturbation A, and the output 
error e are normalized to 1 with all weightings and 
scalings absorbed into the generalized plant struc- 
ture P. This arrangement results in unit invariant 
conditions for robust stability and performance ex- 
pressed in terms of ft that are presented later in this 
section. For the analysis the controller K can be 
considered an (dement of a larger plant G and can 
be absorbed along with the generalized plant P into 
its structure. The diagram for the analysis reduces 
to that in figure 1(b). The analysis problem itself 
involves determining whether e remains in a desired 
set for sets of (/ and A. The resulting structure of G 
c an be partitioned as 


w 


G\\ 

G \2 

z 

V 


G. n 

G22. 

d 


Closing the upper loop of G with the uncertainty 
matrix A results in a linear fractional transformation 
given by 


The conditions for stability and performance ex- 
pressed in terms of H ;3C bounds on portions of the 
generalized plant are given and discussed in some 
detail in references 12 and 13. In the absence of 
uncertainty, the nominal performance objectives are 
expressed in terms of 

IIC22II-X; = sup <r[G 2 2(jw)] < 1 ( 4 ) 

which relates the response c(||c|| 2 < 1) to the set of 
exogenous inputs d(||d|| 2 < 1). In practice, the use of 
scalings and weightings is necessary to represent and 
normalize the varying frequency and spatial content 
of input and output sots. 

Consider plant perturbations that can destabilize 
a nominally stable system. Robust stability is sat- 
isfied for unstructured uncertainty if and only if the 
following condition holds: 

IK'nIlx, < 1 (for all A.ct(A)< 1) (5) 


In general for practical problems, the uncer- 
tainty consists of parameter variations and multiple 
norm-hounded perturbations that result from the im- 
modeled system dynamics. Parameter variations of- 
ten arise from changing flight conditions and rep- 
resent changes in the coefficients of the equations 
that describe the physical system. Unfortunately, 
the norm-bounded test is insufficient and inadequate 
in dealing with robust performance and with realistic 
models of plant uncertainty that involve structure. 

To handle hounded structured uncertainty, the 
structured singular value concept and the function // 
are used to develop necessary and sufficient condi- 
tions. (See ref. 13.) Because any linear problem can 
be rearranged into the form presented in figure 1(a), 
the uncertainty A will be a member of the set defined 
as 

A = |diag(6i , *2 *>m, Ai , A 2 , . . . , A„)|^ € C, 

xAj e C k: J xk 'j | (6) 


( — F ti (G\A)d = \ g *22 + G 2 iA(/ - G 11 A) 1 G 1 2 ] ^ (2) and its bounded subset 


The structure for synthesis is similarly given in fig- 
ure 1(c). The inputs and outputs associated with A 
are absorbed into the exogenous input d! and error 
c'. The equation relating the generalized inputs to 
outputs is given by 

r = f)(P, K)d l = [P n 4- P V2 K[I - P 22 AT 1 Pn]dt (3) 


BA = [A e A|<r(A) < 1] (7) 

where prefix B denotes a unit ball. The function fi 
is defined as 




1 

a I € , AN a )I<^(/-a/a)=oi 


(for A/ e C nxn )(8) 


2 



unless no A e A makes I + M A singular, in which 
case fi(M) = 0. Based on the definition of (i and its 
properties, the robust stability condition for struc- 
tured uncertainty is derived. Robust stability is sat- 
isfied if and only if 

I Km i|| ; , = sup //Km i (jM] < 1 ( for all A e BA) (9) 

Note that, in contrast to <7, the value of fi depends on 
G ii as well as on the structure of the perturbations 
A. In addition, the properties of // permit the 
robust stability of a system with ||Gm iIU = & to 
maintained for all. 

The question of interest to the control designer is 
how well the system performs in the presence of un- 
certainty. Consider the issue of robust performance, 
which describes performance with noise and pertur- 
bations occurring simultaneously. Thus, F U (G. A) is 
stable and 

\\Fn(Cl A)||x. < 1 (for all A € BAifF \\G\\ fl < l) (10) 

The definition of //, is typically not very useful in 
/i computation; however, the extreme cases for the 
structure of A provide the basis for the computa- 
tional bounds on /x. In these special cases in which 
A belongs to one of the extreme sets described below, 
fi is exactly either the spectral radius or the maxi- 
mum singular value 

A = {fi i / n |«€C}-^/i(G)=p(G) 

A = C” X7/ —► //(G) = <r(G) (11) 

Combining these special cases with the fact that // is 
invariant under scaling gives the bounds on ft. Let 
D and Q be the scaling matrices belonging to the set 
defined as given by 

D — {diag(Di , . . . , D rn , d\ /, . . . , d r ,I ) : > 0, dj > 0} 

Q = {Q6*:Q H Q= I Tl } (12) 

Then the bounds on // are 

max p(QG) < //(G) < inf a(DGD (13) 

QgQ De D v / 

The structure of A in general consists of repeated 
scalar and full matrix blocks. The measure // can 
be computed exactly from the upper bound if the 
structure of A corresponds to 25 + F < 3, where 5 
is the number of repeated scalar blocks and F is the 
number of full blocks. Thus, 

/i(G)=inf aiDGD- 1 ) (14) 

Den v ' 


For the problem considered in this paper, A consists 
of a single repeated scalar block and is given in a later 
section. The transformation DGD~ ] is essentially 
a scaling of the inputs and outputs of G and does 
not change the value of //, as mentioned above. In 
addition, because /x can be computed exactly as a 
plus scaling, the methods developed for H ;0C optimal 
control can be used to optimize ft. 

The // analysis can be combined with the H^c 
optimal control to produce /x synthesis to provide 
H^o performance in the presence of structured un- 
certainty. The scaling matrices D and D~ ] are 
used to reflect the structure of A over the frequency 
range. The problem now becomes reformulated as an 
// x -n()nn minimization, 


|| DF,(f\K)D l \\^ < 1 (15) 

known as D K iteration. As the D I\ name implies, 
the ft synthesis approach is to iterate between D 
and K until the solution converges. (See ref. 13.) 
The method is not guaranteed to produce a global 
minimum or to converge; however, the results, widely 
published in the literature, have been successful in 
practical applications. 


Hypersonic Vehicle Model 


A conical accelerator configuration was the ex- 
ample for a generic air-breathing hypersonic vehicle. 
(See ref. 14.) The mathematical model associated 
with the aerodynamic forces and moments, propul- 
sion system, and rigid-body mass moments of iner- 
tia for this vehicle are given in reference 14. An 
updated version of the propulsion model is used in 
this study and is described briefly in appendix A. 
As the vehicle accelerated through a Mach num- 
ber of 8 at an altitude of 85 700 ft, a 10-state lin- 
ear model representing the vehicle dynamics was ob- 
tained at this nonequilibrium flight condition, which 
is characterized by nonzero translational and rota- 
tional accelerations. The conceptual methodology 
for numerical derivation of the linear time- invariant 
models from the six-degree-of- freedom simulation is 
described in appendix B. The linear model was de- 
coupled into a five-state longitudinal and a five-state 
lateral-directional model. The five-state longitudinal 
variable x and control effector u utilized in this study 
arc given by 


n 

a 


x = 


q 

e 

.h. 


(16) 
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The state and control variables are perturba- 
tion quantities and represent deviations from nom- 
inal flight conditions. The nominal flight conditions 
and trim accelerations are provided for reference in 
table I. The open- loop characteristics of the plant are 
unstable. The eigenvalues arc as follows: 

Open- Loop Eigenvalues 


at hypersonic speeds with a change 1 in angle of at- 
tack. The original nonlinear model assumes that rh a 
is an instantaneously adjustable quantity and thus is 
kept constant with changing flight conditions, which 
results in %Aij = %A rhj. Therefore, to reflect real- 
istic utilization and to forgo time-consuming changes 
in the nonlinear model, the control problem is formu- 
lated with the fuel flow rate as the control effector 
and its changes are presented as the percentage of 
change from the nominal rather than as actual units 
of mass per second. 


<>.5-103 x 10 

- 9.13-1 x K) 1 ± j3.557(i x 10" 2 
-1.8312 
1 .6533 

Note that the altitude is included as a state variable 
to acc ount for temperature, density, and gravity gra- 
dients. These variations significantly affect the longi- 
tudinal long-period dynamics of the vehicle and add 
an aperiodic altitude mode caused by the variation 
of atmospheric density with altitude. (See ref. 15.) 

Several interesting nuance's about this model are 
worth considering. In all practical problems, a lin- 
ear time invariant. (LTT) system is only an approxi- 
mation for the real physical system. In most cases, 
time invariance 1 of system characteristics is valid for 
small increments of time dictated by the problem. In 
this particular case, the vehicle accelerates through 
a Mach number of 8 and the LTI system is valid for 
only that instant. However, if the parameter varia- 
tion with time is represented as a multiplicative un- 
certainty such that G ~ G(I + A), then the linear 
system can be considered time invariant over some 
interval and can be used in LTI centred design. In 
this initial application of the and // control the- 
ories, all uncertainty is assumed to be represented 
by indivielual uncertainty on the effectiveness of the 
eleven and the filed flow rate. This representation 
fulfills two desirable guideline's it keeps the initial 
application uncomplicated and it introduces struc- 
tured uncertainty into the problem. 

The second important point is fuel equivalence 
ratio as a control effector. The fuel equivalence ratio 
if is defined as 


where ihj is the fuel mass flow rate, / is the stoi- 
chiometric ratio that for hydrogen is 0.029, and rh a 
is the air mass flew rate. In practice, it is important 
te) recognize that m n changes almost instantaneously 


The' dependence of inlet conditions on angle-of- 
attack variation is not explicitly included in the non- 
linear model. However, the control problem is formu- 
lated to consider the major aerodyuamic/propulsive 
interactions. Thc?sc interactions are reflect cel in ex- 
plicit performance requirements on angle of attack 
as well as on the engine control effector bi'uj* rather 
than on the fuel (equivalence ratio. Using // as a con- 
trol effector implies that either the angle' of attack 
has no effect on m n flow into the propulsion system 
or that m a is instantaneously measurable, thus allow- 
ing for matching adjustment to ihf to get the desired 
?/. However, neither implication is realistic. Hence, in 
practice, while' the system is following a commanded 
//, tii'nj must be used as a control effector. 

Robust Control Problem Formulation 

The controller requirements were established 
based on the near-optimal ascent trajectory and on 
the sensitivity of the air-breathing propulsion sys- 
tem to angle-of-attack variations. Thus, the speci- 
fications in the presence of atmospheric turbulence 
and system uncertainty are the following: highly ac- 
curate tracking of velocity and altitude commands, 
limiting angle-of-attack total deviation to «D.f>° from 
nominal, minimizing control power use, and stabi- 
lizing the vehicle. For this study. x sVti \v = Vinous 
was assumed available for output feedback; the as- 
sumption was based on the availability of inertial 
angle-of-attack calculations with feedback from an 
inertial measurement, unit (IMU). The reason for in- 
ertial angle of attack is twofold. First, a conventional 
angle-of-attack measurement that uses a pitot tube 
in a free stream is not practical at hypersonic speeds 
because of very high stagnation temperatures. Sec- 
ond, angle-of-attack perturbations from atmospheric 
turbulence are significant compared with the desired 
response (^20 percent); a signal of such magnitude 
should not be introduced directly into the feedback 
loop. Measuring the aerodynamic angle of attack 
would introduce atmospheric turbulence into the con- 
trol problem as measurement noise rather than as a 
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disturbance, thus changing the way it would be han- 
dled by the control law. 

The block diagram problem formulation is illus- 
trated in figure 2(a). All feedback state measure- 
ments were assumed corrupted by noise. The noise 
matrix UV.oisr- was represented by 

ww = 10 ( 19 ) 

The noise matrix was only partly intended to rep- 
resent realistic sensor data; rather, it was included 
to satisfy current computational requirements associ- 
ated with solving the //-*, problem. The control actu- 
ator dynamics were represented by first-order filters 
with a 30 rad/sec bandwidth for the elevon and a 100 
rad/sec bandwidth for the fuel flow rate. In a generic 
air-breathing hypersonic vehicle, the uncertainty, as 
discussed in the “Introduction^ comes from different 
source's and occurs simultaneously. Thus, the physics 
of the problem imposes a structure on the uncertainty 
of a hypersonic vehicle. For this problem, system un- 
certainty due to parameter variation arising from ve- 
hicle acceleration was represented as variations in the 
effectiveness of the commanded elevon and the fuel 
flow rate. Specifically, a 20-pcreent multiplicative 
uncertainty in control effectiveness was introduced 
where 


conical accelerator. Thus, frequency dependent per- 
formance weightings are used for velocity and al- 
titude errors. Such time history requirements as 
percentage of overshoot , steady-state error, and re- 
sponse time constant are translated into first-order 
transfer functions given by equations (22) and (23). 
The details for deriving these transfer functions as 
well as for other performance weightings are given in 
appendix C. Thus, the velocity performance is spec- 
ified by 


0.5(.s + 4.330 x 10 2 ) 
,s + 2.165 x 10 - V- 


( 22 ) 


Similarly, altitude response' is constrained by 




O.T)(,s -i- 3.819 x l(r 2 ) 

.sTV.b^'Tx ’ll) 1 


(23) 


The performance weightings on velocity and alti- 
tude' describe; the response of each variable to the 
commanded value and to the' ability to reject- dis- 
turbance's. Low-frequency gain elescribe*s both the 
steady-state error and the magnitude of disturbance 
rejection, which is discussed in more detail in appen- 
dix C. Because velocity and altitude were the only 
variables tracked, constant weighting was applied to 
angle of attack a, e/, and 0. resulting in 


A = 


h 

0 


0 

h 


(#€C||A||^< 1) (20) 


and the uncertainty weighting matrix was 


W A = 


0.20 

0 


0 

0.20 


(21) 


Thus, the control effectiveness was forced to vary 
from 80 to 120 percent of the assumed nominal value. 
At this stage of problem development, atmospheric 
density perturbations were assumed reflected in the 
20-percent uncertainty in fuel flow rate effectiveness. 

Performance specifications for a flight control sys- 
tem translate quite well into an Hqq context for this 
problem. In designing for tight performance margins, 
performance weighting functions were augmented to 
the system. Performance specifications for Hqc make 
practical sense only when meaningful variables are 
specified for weighting functions. As illustrated in 
figure 2(a), the weighting functions for control effec- 
tor positions and rates as well as for the state vector 
were employed as output performance variables. 

Specific performance requirements were derived 
from a near-optimal fuel ascent trajectory for the 


a 


"20' 

<} 

= 

1 

_ 0 _ 


1 


(24) 


The weighting for a was based on the desire to at- 
tenuate atmospheric disturbances and to limit the 
magnitude of the output as much as possible with- 
out violating the performance requirements of other 
variables. The responses of q and 6 , which are in- 
directly related to n response, were not of primary 
interest, as reflected in the weighting selection. The 
actuator position and rate limits were imposed by 


Wp 


he 


30 

hihf 


60 

he 


20' 

>»/. 


10 


(25) 

(26) 


Atmospheric turbulence is also explicitly included 
in problem formulation. The primary concern is the 
effect of turbulence on engine performance. Turbu- 
lence can affect performance directly either by chang- 
ing inlet flow conditions or by exciting actuator con- 
trollers; in both cases, the result will be undesirable 
vehicle angle-of-attack variations. Longitudinal and 
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vertical Dryden turbulence filters are implemented as 
input weighting functions in the generalized frame- 
work (F u and F lv in figure 2(a)). (See refs. 16 and 17.) 
Filter inputs w \ and iry to F u and F w . respectively, 
represent white zero-mean unit variance noise sig- 
nals. The longitudinal Dryden filter output is as- 
sumed to act along the velocity vector due to a small 
angle-of-attack flight condition. The vertical Dry- 
den filter is divided by nominal velocity V 0 to give 
the angle of attack instead of the vertical velocity 
perturbation. 

The block diagram in figure 2(a) can be manip- 
ulated into the general framework of figure 1(a) as 
depicted in figure 2(b). Recall that all the input 
and output signals of the generalized plant P. ex- 
cept for the control input and measured output, be- 
long to the unity bounded sets with scaling absorbed 
into P. In this problem the performance weight- 
ing functions also served as the scaling factors for 
the output signal set. Thus, the input labels in fig- 
ure 2(b) refer to the physical quantities represented 
by the input and output in figure 2(a). Note that 
the controller-commanded input data are also, after 
weighting, the matrix input data and thus define 
the system uncertainty. For control synthesis, the 
commands, noise, and atmospheric turbulence (con- 
tained in d) and the actuator uncertainty (contained 
in z) are combined into augmented d as denoted in 
figure 1(c). The results of the controller design pro- 
cess are discussed in the next section and the nu- 
merical results for the state-space problem formula- 
tion and controller poles and zeros are provided in 
appendix D. 

Robust Control System Design Results 
and Evaluation 

Three different controllers are designed and eval- 
uated for the problem formulated in the prior section. 
These controllers include the controller designed 
lor a nominal system that contains only measurement 
noise to corrupt the feedback signal. The second de- 
sign is again the H x controller designed explicitly 
to tolerate 20-percent uncertainty in control effec- 
tiveness. The final design is a ft controller based on 
the second // x controller. The designs are individ- 
ually evaluated and compared for how successfully 
they provide the desired level of performance in the 
presence of the prescribed level of uncertainty. The 
dynamic characteristics of the controllers are given 
in appendix D. 

As the initial step in a controller design, an //^ 
controller for a model with no uncertainty is ob- 
tained. The nominal airplane model is derived from 


the system interconnection shown in figure 2(b) ei- 
ther by deleting the rows and columns of P corre- 
sponding to w and 2 or by setting - 0. Fre- 
quency domain closed-loop system analyses of this 
// x controller for nominal performance, robust sta- 
bility. and robust performance are illustrated in fig- 
ure 3. Robust stability and robust performance wort' 
measured for the structured uncertainty using //.(•) 
tests of the appropriate transfer function matrix. 
(Sec the section, “Theoretical Review of //.") The 
closed-loop system that uses the // x controller sat- 
isfies nominal performance requirements for simul- 
taneous inputs of two tracking commands, velocity 
and altitude, and in the presence of longitudinal and 
vertical atmospheric turbulence. Recall from equa- 
tion (4) that nominal performance is satisfied if and 
only if o\G'>2{ju)] < 1 for all frequencies illustrated 
in figure 3. Satisfying the nominal performance con- 
dition signifies that specified response characteris- 
tics are met for the worse possible combination of 
bounded inputs into the nominal system (i.e., no un- 
certainty). The lack of robust stability and robust 
performance in the presence of uncertainties is not 
surprising here because the uncertainties wore not 
explicitly taken into account. This // x controller is 
used as the baseline in the subsequent, comparisons 
of controller performance. 

I he time histories of the closed-loop system re- 
sponse were obtained from a linear simulation in 
this initial application of // x and \i robust con- 
trol theories. The linear time simulation was used 
to translate 1 frequency domain results into the time 
domain, which is more conducive to evaluating ac- 
tual physical system performance. As mentioned 
in the problem formulation, moderate atmospheric 
turbulence was implemented through lough udinal 
and vertical Dryden filters. Moreover, for consis- 
tent implementation in a discrete simulation envi- 
ronment, the continuous-time Dryden filters were 
divided by \J Ny c \ 1 1 i s t frequency . The atmospheric 
turbulence implementation is discussed in more de- 
tail in appendix E. 

The time response of the nominal airplane model 
to simultaneous commands of a velocity change of 
100 ft/sec and an altitude change of 1000 ft in the 
presence of longitudinal and vertical atmospheric tur- 
bulence is presented in figures 4(a) 4(e). These fig- 
ures also include the fi controller nominal perfor- 
mance that will be discussed later in this section. 
Both velocity and altitude (figs. 4(a) and (b)), meet 
the performance requirements derived from track- 
ing a near-optimal fuel trajectory. It is interesting 
to note that the absence of uncertainty in the con- 
troller design did not benefit nominal performance 
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as may have been suspected. The lack of uncertainty 
in eleven effectiveness made this control surface too 
sensitive to atmospheric turbulence as reflected in the 
second-order behavior superimposed on the primary 
response of the altitude as the commanded value is 
reached. Even when the altitude variable displays 
this response, the nominal performance' criterion is 
satisfied. Changing the performance weighting on ei- 
ther the eleven or the altitude does not satisfactorily 
resolve this problem. The propulsion system effective 
angle of attack shown in figure 4(c) also fulfills per- 
formance' spec ifications for total deviation of ~0.T)°. 
The eleven total deflection in figure 4(d) is roughly 
1.5°. thus limiting control-surface-induced drag. The 
fuel flow rate' shown in figure 4(e) avoids large sudden 
changes in magnitude', t hus minimizing transients in 
the combustor. 

If cont rol surface deflection as well as the magni- 
tude' and rate of fuel flow rate can be reduced, ve- 
hicle' performance will improve. Typically, for an 
air-breathing SSTO vehicle, the payload fraction is 
«3 percent, and the fuel fraction is « (>0 percent. 
(Sec* ref. A.) Hence, any improvement in fuel frac- 
tion due to reduction in control-surface-induced drag 
has a potential impact on the consumed fuel that is a 
large percentage of vehicle weight. Furthermore, the 
performance of the propulsion system is extremely 
sensitive' to (“hanging conditions in the inlet, com- 
bustor, and nozzle; therefore, transients in all parts 
of the engine should be minimized. The control sys- 
tem minimizes perturbations in the inlet conditions 
by limit ing angle' of attack and aids combustion sta- 
bility with smooth change's in fuel flow rate. 

As mentioiH'd, the design goal is to maintain per- 
formance with a spe'cified 20-percent control power 
uncertainty. The baseline 12-state H oc. controller, 
discussed above, is used for c losed-loop system anal- 
ysis. Before robust performance can be considered, 
robust stability must be addressed. As is evident 
from figure 3, the robust stability requirement is not 
satisfied by the Hoc controller for the specified uncer- 
tainty. Recall from the section “Theoretical Review 
of //’* that the dosed- loop system is internally stable 
for all perturbations with magnitude l//i where 8 is 
the peak magnitude of the frequency response. For 
a system with inputs normalized to 1, as is done in 
the Hoc control problem, A — 1 for 1/8. If 8 > 1 , 
the implication is that the system is destabilized by a 
perturbation with a smaller magnitude than the spec- 
ified uncertainty. The fact that fi[G\ \ (ju;)] — 1.174 
implies that the closed-loop system remains stable 
only for 17-percent uncertainty in the control effec- 
tiveness of both the elevon and the fuel flow rate. 
The details are elaborated upon in Appendix D. The 


robust performance is also not satisfied by the H oc 
controller because it requires robust stability as a 
necessary condition. 

An H^, controller that is designed with 20- 
perccnt uncertainty explicitly included in the system 
(it is treated as unstructured by Ho c optimization) 
does not fulfill performance requirements even for an 
ideal system. In fact, the nominal performance con- 
dition ||G22 ||oc < 1 violated as illustrated in fig- 
ure 5. Figure' (i provides a sample time response' of 
velocity to a 100 ft /sec step command. The time rc'- 
sponse clearly shows that the* velocity response does 
not achieve the' desired steady-stale' error (<5 per- 
cent) or the' rise' time' (40 se'c). In fact, the response 
to the velocity command violates the nominal per- 
formance criterion. Further analysis indicates that 
either system uncertainty conditions or performance? 
specifications 011 the* tracking variables must be ro- 
laxe'd; otherwise, an Hoc controller cannot, fulfill ro- 
bust performance requirements as specified for this 
problem. However, this cemtroller serves as a basis 
for computing a // controller. 

In an attempt to satisfy both robust stability 
and robust performance with original specifications, 
a // controller is computed based on D K iterations. 
The resulting /e controller is reduced from 18 to 13 
states by high-frequency rosidualized truncation and 
Hankel optimal norms. The nominal performance of 
the // controller compares we'll with the nominal per- 
formance? of the* baseline Hoc controller as illustrated 
in figures 4(a) 4(e). Note' the lack of the superim- 
poses l second-order response 011 the' altitude' variable' 
for the /1 controller. However, in the robust, per- 
formance, the advantage' of a // controller becomes 
apparent. 

The frequency domain analysis indicates that the 
performance requirements are satisfic'd in the pres- 
ence of 20-percent control effector uncertainty, so 
that the robust performance' condition in figure 7 is 
met. The velocity, altitude, and angl('-of-attack time 
histories are essentially the same for nominal and per- 
turbed systems, as illustrate'd in figuies 8(a) 8(e). 
The noticeable deviation of the perturbed response 
from the nominal value occurs for the variable's that 
contain uncertainty. The time histories of the con- 
trol variable's exhibit differences that are due to the 
positive or negative value of the uncertainty because' 
the controller tries to compensate for the uncertainty 
in effectiveness while attempting to achieve the de- 
sired response. The elevon response in figure 8(d) 
shows a small variation in the magnitude of the ini- 
tial deflection; a smaller than nominal deflection for a 
20-pcrcent increase in effectiveness; and, conversely, 
a larger deflection for a 20-percent deterioration of 
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effectiveness. The fuel How rate in figure 8(e), shows 
that the difference between nominal and perturbed 
responses is more pronounced, although smoothness 
still characterizes each response. 

In fact, the robust performance margin suggests 
t hat the system could withstand higher levels of con- 
trol effector uncertainty and still fulfill performance 
requirements. The frequency analysis shown in fig- 
ure 9 indicates that robust- performance can be sus- 
tained for a closed-loop system with a 32. 5-percent 
control effector uncertainty. The A used in the? time 
simulations is a worst case, real- rational, stable per- 
turbation with IIApllx, — 1.0. The closed-loop sys- 
tem time response is provided in figures 10(a) 10(e) 
for the nominal system and the ±32. 5-percent uncer- 
tainty in control effectiveness. A comparison of nom- 
inal and worst case perturbed time responses shows 
that the // controller successfully handles actuator 
unc ertainty without sacrificing performance in most 
system variables. Some visible performance degra- 
dation occurs in the variables directly affected by 
the uncertainty, eleven, and fued flow rate, but this 
degradation is not surprising. Despite some perfor- 
mance degradation compared with the ideal system 
response, all requirements on every system variable 
are satisfied, even with higher uncertainty than the 
20 percent proscribed in the design. The encourag- 
ing preliminary results of the // controller establish 
the technique as potentially successful in deviling with 
unique characteristics of hypersonic vehicles. 

Conclusions 

The applicability of robust control techniques to 
a single-stage- to-orbit air-breathing vehicle at hyper- 
sonic speeds on an ascent accelerating path and the 
effectiveness of these techniques are explored in this 
paper. Uncertainty is an integral part of hypersonic 
vehicle characteristics, so its effect on the analysis 
and synthesis of various control system design tech- 
niques is important to understand. 

Several important issues related to control sys- 
tem design should be noted. The characteristics 
of air-breathing hypersonic vehicles and the require- 
ments imposed on the control system translate ex- 
plicitly into domain specifications. Atmospheric 
turbulence has a noteworthy effect on the design of 
the //x: controller. The best performance should be 
exhibited bv the controller designed for a nominal 
system (i.e., the assumption is that the physical sys- 
tem is perfectly represented by the model). However, 
the very fact that the system is apparently perfectly 
known makes the elevon sensitive enough to atmo- 
spheric turbulence to induce residual second-order 
behavior in the altitude response that it directly con- 
trols. In addition, this nominal H 0 Q controller does 
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not remain stable at 20- percent uncertainty, which is 
the prescribed level for this design problem. 

The //x controller, designed explicitly for 
20-percent uncertainty, failed to provide the desired 
level of nominal system performance. The // con- 
troller, designed for the same level of uncertainty, 
preserved the required performance and tolerated 
more than 50 percent more uncertainty than the de- 
sign specified, thereby providing robust per fori nance. 
This initial application indicates that the nature of 
uncertainty and how it, is handled by the control 
design methodology have significant effects on the 
achievable level of system performance as well as on 
the tolerable level of system uncertainty. 

Frequency- based linear analysis technique's imply 
the importance of fi as both the analysis and the 
synthesis tool for an air-breathing hypersonic vehi- 
cle. Even though the control laws were not tested 
in a full six-degree-of- freedom nonlinear simulation, 
the results are promising after the initial applica- 
tion of the structured singular value theory to an 
air-breathing hypersonic vehicle. The relatively fast 
design time is conducive to the tradeoff study be- 
tween the achievable level of performance and the 
prescribed level of uncertainty. Such a study is valu- 
able for the conceptual design as well as for the actual 
vehicle because the known physical parameter varia- 
tions are explicitly considered and their relationships 
to the physical system are preserved. A timely trade- 
off study can determine how much of the true phys- 
ical uncertainty needs to be directly represented in 
the controller design to achieve the desired level of 
performance for the prescribed level of uncertainty. 

Furthermore, because uncertainty occurs simulta- 
neously from many different sources and the degree of 
uncertainty is high, the physical characteristics of the 
air-breathing hypersonic vehicle introduce structure 
into the problem. A methodology that takes full ad- 
vantage of these physical characteristics is essential. 
The //-analysis and synthesis technique preserves the 
structural relationship between uncertainty and per- 
formance variables, allowing the designer a system- 
atic approach to explore tradeoffs between the two. 
Failure to account for this structural relationship can 
result in excessively conservative specifications and 
poor designs for an air-breathing hypersonic vehicle. 
Although further study is necessary before the struc- 
tured singular value technique can be recommended 
as the method of choice for air-breathing hypersonic 
vehicles, the beginning research has been promising. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
April 1, 1994 



Table I. Nominal Flight Conditions and Trim Accelerations 

Mach 8*0 

Dynamic pressure, psf 2005.8 

Weight, lb 261 125.9 

Moment of inertia, I yy , slug- ft 2 9202.797 

Center of gravity, x cg , ft 10.0773 

V , ft /sec 2 30.6172 

а, deg/sec —1.4309 x 10 

< 7 , deg/sec 2 5.3973 x 10 

б , deg/sec 1.634 x 10 2 

A, ft/sec 163.146 

V , ft/sec 7851.6 

a, deg 2.1628 

g, deg/sec 3.7754 x 10 

0, deg 3 3534 

M t 85 745.7 

Elevon Se Q , deg —9.0720 

Throttle setting, r] 0 1.2743 

Fuel flow rate, lbm/sec 171.45 
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(a) General structure. 



(b) Analysis structure. 



(c) Synthesis structure. 

Figure 1. General framework used in robust, feedback control system design. Any linear combination of inputs, 
outputs, and commands along with perturbations and controller can be viewed in context and rearranged 
to match appropriate diagram. 
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(b) Physical description of unity-bounded input and output sets of P. 
Figure 2. Conical accelerator example problem formulation. 
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Figure 3. Frequency analysis of closed-loop system with H ^ controller perturbed by 20-percent uncertainty in 
control effectiveness of both elevon and fuel flow rate. 
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(a) Velocity. 



(b) Altitude. 


Figure 4. Time response of nominal closed-loop system commanded to simultaneous 100 ft /sec ve 
1000-ft altitude change while subjected to vertical and longitudinal atmospheric turbulence. 
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(c) Angle of attack. 











Figure 5. Nominal performance frequency analysis of closed-loop system with controller designed for 
20-percent uncertainty in control effectiveness of both elevon and fuel flow rate. 
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Figure 6. Velocity response to velocity command, altitude command, and vertical and longitudinal turbulence 
of system with H 0 c controller designed explicitly with 20-percent uncertainty in control effectiveness of both 
elevon and fuel flow rate. 
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Figure 7. Frequency analysis of closed-loop system with /x controller perturbed by 20-percent uncertainty in 
control effectiveness of both elevon and fuel flow rate. 
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(a) Velocity. 



(b) Altitude. 


Figure 8. Time history of nominal and perturbed (±20-percent uncertainty in effectiveness of elevon and fuel 
flow rate) closed-loop system with (i controller commanded to simultaneously increase velocity by 100 ft/sec 
and altitude by 1000 ft while subjected to vertical and longitudinal atmospheric turbulence. 
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ure 9. Frequency analysis of closed-loop system with /x controller perturbed by 32.5-percent uncertainty 
control effectiveness of both olevon and fuel flow rate. 




(a) Velocity. 



(b) Altitude. 


Figure 10. Time history of nominal and perturbed (±32. 5- percent uncertainty in effectiveness of elevon and fuel 
How rate) closed-loop system with ft controller commanded to simultaneously increase velocity by 100 ft/sec 
and altitude by 1000 ft while subjected to vertical and longitudinal atmospheric turbulence. 
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Appendix A 

Upgraded Propulsion Model for Conical 
Accelerator 


for a specified Mi and 771 . Equation (Al) can be 
derived using the definition of /sp and r] and the as- 
sumption that the air mass flow through the engines 
is a constant at a given Mach number. Thus, 


The upgraded propulsion model for a conical ac- 
celerator presented in this appendix replaces the one 
given in reference 14. It was generated to eliminate 
certain inconsistencies in the original database and to 
provide a gridded, smooth propulsion model suitable 
for trajectory simulation, analysis, and optimization. 
The propulsion model in reference 14 had inconsis- 
tent characteristics in estimating propulsion system 
performance at values of rj away from unity. These 
inconsistencies occurred because the propulsion mod- 
eling codes were run at relatively few analysis points 
(about 7 / — 1 for different Mach numbers) and t he 
database was then padded with extrapolated values 
to extend the range. However, the extrapolated data 
points violated a fundamental relationship between 
thrust coefficient (Cj\ ft 2 thrust per dynamic pres- 
sure) and specific impulse. 

The upgraded model makes use of the data at 
certain calculated points that were theoretically pre- 
dicted (ref. 14) and approximates the generic per- 
formance of this class of propulsion systems at grid 
points where the data are unavailable. At a given 
Mach number, Cj was assumed to have a nearly lin- 
ear variation with ?/ for values of ?/ exceeding unity. 
Moreover, Cy was assumed to rapidly decrease as 
// was decreased below 0.3 to simulate an engine 
unstart. A general curve- fitting algorithm was em- 
ployed to fit Cj data at certain calculated points 
while maintaining the aforementioned trends. The 
curve fit was then used to extract Cy at a prespecified 
set of 7 / values. To obtain the specific impulse (^spi 
see) at this set. of // values, the following relationship 
was used: 


/s P (r/.Mi) = 


Criv, Mi)isy{vu M\)m 
CrtiuMi) 


(Al) 


7h f 

0.029 rh a 


(A2) 


/sp = ^ (A3) 

m f9 

In the above equations, q and <j are dynamic pressure? 
and acceleration due to gravity, respec tively. The 
above procedure was repeated at all Mach numbers to 
generate’ a grid. Furthermore, tlu' upgraded propul- 
sion inodefl is simplified by eliminating the variation 
of C T and / s p with q and by taking into account the 
weak dependence, especially in the lower Mach num- 
ber range (A/ < 15). 

The upgraded model consists of two sets of tables: 
one set for the low-speed propulsion cycle (Mach = 
0.3 to 2) and another set for the high-speed propul- 
sion cycle (Mach = 2 to 25). The Gy and / S p val- 
ues are given in tables Al and All for the low-speed 
propulsion cycle and in tables A I II and A IV for the 
high-speed propulsion cycle. Axial thrust T (lb) and 
fuel flow rate lT(lb/sec) are computed as follows: 

T = qCy(hLq) (A4) 


W 


T 

MA/.//) 


(A5) 


Note that the propulsion model is not continu- 
ous at a Mach number of 2, which would repre- 
sent discrete switching from low- to high-speed cy- 
cles. To avoid this discontinuity, the data could 
he linearly interpolated between a Mach number of 
1.5 (low-speed cycle) and a Mach number of 2.0 
(high-speed cycle) to simulate a gradual transition 
from the low- to the high-speed propulsion cycle. 
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Table AI. Low-Speed Thrust Ooeffieent as Function 
of Mach Number and Fuel Equivalence Ratio 


Low-speed thrust coefficient for Mach number of 


n 

0.3 

0.5 

0.7 

0.9 

1.0 

1.5 

“ 2.0 ' 

0.5 

1047.6 

503.61 

364.83 

294.02 

273.68 

198.56 

163.28 

.6 

1257.0 

595.11 

425.86 

341.86 

318.28 

233.55 

194.01 

.7 

1466.5 

683.53 

482.91 

386.03 

359.50 

266.97 

224.09 

.8 

1675.9 

768.88 

535.98 

426.54 

397.34 

298.82 

253.54 

.9 

1885.3 

851.15 

585.08 

463.39 

431.81 

329.10 

282.34 

1.0 

2094.7 

930.35 

630.19 

496.58 

462.89 

357.81 

310.50 

LI 

2301.0 

1006.5 

671.32 

526. 1 1 

490.59 

384.94 

338.02 

1.2 

2513.3 

1079.5 

708.48 

551.98 

514.92 

410.50 

364.89 

1.3 

2722.6 

1149.5 

74 1 .65 

574.19 

535.87 

434.49 | 

391.12 

1.4 

2931.9 

1216.4 

770.85 

592.75 

553.44 

456.91 

116.72 

1.5 

3141.2 

1280.2 

796.06 

607.64 

567.63 

477.75 

411.66 

1.6 

3350.5 

1341.0 

817.30 

618.88 

578.44 

497.02 

465.97 

1.7 

3559.7 

1398.6 

834.56 

626.45 

585.87 

514.72 

189.63 

1.8 

3768.9 

1453.2 

847.84 

630.37 

589.92 

530.85 

512.65 

1.9 

3978.1 

1504.7 

857.14 

630.63 

590.58 

545.40 

535.03 

2.0 

4187.2 

1553.2 

862.46 ; 

627.23 

587.89 

558.38 

556.77 


Table AIL Low-Speed Specific Impulse as Function 
of Mach Number and Fuel Equivalence Rat io 


Low-speed specific impulse for Mach number of 


n 

0.3 

0.5 

' 0.7 T 

0.5 

1894.0 

2332.2 

2896.3 

.6 

1893.9 

2296.6 

2817.3 

.7 

1893.8 

2261.0 

2738.3 

.8 

1893.7 

2225.4 

2659.4 

.9 

1893.6 

2189.8 

2580.4 

1.0 

1893.5 

2154.2 

2501.4 

1.1 

1893,4 

2118.6 

2422.5 

1.2 

1893.4 

2083.0 

2343.5 

1.3 

1893.3 

2047.4 

2264.5 

1.4 

1893.2 

2011.8 

2185.5 

1.5 

1893.1 

1976.2 

2106.6 

1.6 

1893.0 

1940.6 

2027.6 

1.7 

1892.9 

1905.0 

1948.6 

1.8 

1892.8 

1869.4 

1869.6 

1.9 

1892.7 

1833.8 

1790.7 

2.0 

1892.6 

1798.2 

1711.7 


0.9 

1.0 

1.5 

2.0 

3254.9 

3299.5 

2860,1 

2512.6 

3153.7 

3197.7 

2803.8 

2487.8 

3052.4 

3095.8 

2747.2 

2463.1 

2951.1 

2994.0 

2690.6 

2438.4 

2849.9 

2892.2 

2633.9 

2413.7 

2748.6 

2790.3 

2577.3 

2389.0 

2647.3 

2688.5 

2520.7 

2364.3 

2546.0 

2586.6 

2464.0 

2339.6 

2444.8 

2484.8 

2407.4 

2314.9 

2343.5 

2383.0 

2350.8 

2290.2 

2242.2 

2281.1 

2294.2 

2265.4 

2111.0 

2179.3 

2237.5 

2240.7 

2039.7 

2077.4 

2180.9 

2216.0 

1938.4 

1975.6 

2124.3 

2191.3 

1837.1 

1873.7 

2067.7 

2166.6 

1735.9 

1771.9 

2011.0 

2141.9 
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Table AIII. High-Speed Thrust Coefficent as Function of 
Mach Number and Fuel Equivalence Ratio 




High-speed thrust coe 

fficicnt for Mach number of 



V 

2 

3 

4 

6 

8 

10 

15 

20 

25 

0.25 

205.04 

147.41 

97.613 

5.631 

65.041 

49.924 

28.204 

24.548 

12.904 

.50 

315.57 

236.30 

181.66 

109.21 

113.27 

83.254 

50.136 

44.460 

24.635 

.75 

391.29 

303.64 

252.14 

160.75 

140.52 

102.74 

64.496 

58.078 

33.843 

1.0 

451.44 

361.34 

309.05 

210.23 

156.16 

115.38 

74.624 

68.082 

41.566 

1.5 

553.90 

465.68 

409.85 

307.23 

171.55 

130.85 

88.167 

82.205 

54.488 

2.0 

651.48 

566.99 

510.65 

404.23 

180.42 

140.50 

97.357 

! 92.411 

65.602 

2.5 

753.04 

670.76 

611.45 

501.23 

189.71 

147.95 

104.77 

100.97 

75.82 

3.0 

861.74 

778.94 

712.25 

598.23 

200.98 

154.73 

111.59 

108.96 

85.59 

4.0 

1104.1 

1010.7 

913.84 

792.23 

225.90 

168.77 

125.64 

125.28 

104.68 

5.0 

1379.6 

1263.1 

1115.4 

986.23 

242.09 

184.95 

141.73 

143.51 

123.85 


Table AIV. High-Speed Specific Impulse as Function 
of Mach Number and Fuel Equivalence Ratio 



High-speed specific impulse for Mach number of 

V 

2 

3 

4 

6 

8 

10 

15 

20 

25 

0.25 

5672.4 

5054.6 

3877.4 

2891.8 

3760.2 

3214.1 

1948.7 

1885.1 

1096.5 

.50 

4365.0 

4051.3 

3607.9 

2838.5 

3274.3 

2679.9 

1732.1 

1707.0 

1046.7 

.75 

3608.3 

3470.6 

3338.5 

2785.3 

2708.0 

2204.6 

1485.4 

1486.6 

958.60 

1.0 

3122.2 

3097.5 

3069.0 

2732.0 

2257.0 

1857.0 

1289.0 

1307.0 

883.00 

1.5 

2553.9 

2661.3 

2713.3 

2661.7 

1652.9 

1404.0 

1015.3 

1052.1 

771.68 

2.0 

2252.9 

2430.2 

2535.5 

2626.5 

1303.8 

1130.6 

840.85 

887.02 

696.81 

2.5 

2083.2 

2300.0 

2428.8 

2605.5 

1096.8 

952.45 

723.89 

775.34 

644.28 

3.0 

1986.6 

2225.8 

2357.6 

2591.4 

968.24 

830.11 

642.53 

697.28 

606.11 

4.0 

1909.0 

2166.1 

2268.7 

2573.8 

816.24 

679.04 

542.54 

601.27 

555.95 

5.0 

1908.3 

2165.5 

2215.4 

2563.3 

699.78 

595.35 

489.64 

551.01 

526.19 
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Appendix B 

Linear Model Development From 
Nonlinear Simulation 

This appendix describes the conceptual develop- 
ment of the linear models of a conical accelerator 
vehicle. The linear models were derived numeri- 
cally from the 6-DOF nonlinear rigid-body simula- 
tion ot the vehicle using the program to optimize 
simulated trajectories (POST). (See ref. 18.) The 
POST simulation is for a general, rigid-body, 6-DOF 
trajectory with discrete parameter trajectory target- 
ing and optimization. The 6-DOF simulation of the 
vehicle is built by specifying the vehicle aerodynam- 
ics. propulsion, and mass properties and the environ- 
ment in which the vehicle operates. For the 6-DOF 
rigid-body simulation of a conical accelerator, the 
aerodynamics and mass properties in reference 14 
were used along with the propulsion model given in 
appendix A. To derive the linear models, a spherical 
nonrotating planet, model with a stationary atmo- 
sphere (a 76 U.S. Standard Atmosphere option cho- 
sen in POST) was assumed to define the operating 
environment. Linear models of the vehicle dynamics 
at the desired High! conditions were extracted using 
the targeting feature of the POST simulation. 

The targeting option in the POST simulation al- 
lows a user to solve' a nonlinear programming prob- 
lem wherein a user-defined set of dependent variables 
(constraints) must be driven to their desired value's 
by changing a set of user-defined independent vari- 
able's (controls). The solution of the nonlinear pro- 
gramming problem is obtained from a gradient-based 


optimization algorithm. A special case of the tar- 
geting problem is the trim problem wherein a set 
of dependent variables, usually ve'hicle accelerations 
at a given instance are inquired to take on some 
user-specified set of desired values by allowing the 
program to change' a set of independent variables 
(i.e., control effectors). The trim problem and the 
fact that the solution of this problem is obtained 
bv a gradient- based algorithm are exploited to nu- 
merically derive the linear models at a given flight 
condition. 

To get the linear models, a special trim problem 
is posed in the POST simulation. In this case, the 
dependent variables are defined to be vehicle transla- 
tional and rotational accelerations; the independent 
variables are vehicle states and control effectors. To 
solve this problem, the POST simulation calculates 
the numerical partial derivative of each of the de- 
pendent variables with respect to each independent 
variable by using a forward differencing scheme. The 
Jacobian, which consists of the first-order terms of 
the Taylor series expansion of the nonlinear equa- 
tions of motion about the analysis point (the cho- 
sen flight conditions), is extracted and partitioned 
appropriately to give the linear system state space 
representation. Note that the analysis point is com- 
pletely arbitrary; that is, the analysis point could be 
a nonequilibrium flight condition wherein vehicle ac- 
celerations are nonzero or it could be a steady-state 
flight condition. However, nonlinear equations of 
motion linearized at nonequilibrium flight conditions 
result in linear systems that are time varying and spe- 
cial care must be exercised while using them in eon- 
junction with LTI design and analysis procedures. 
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Appendix C 

Weighting Function Derivation 

The performance weighting function W p enters 
into the control law design through the following 
relationship: 


The desired percentage of overshoot is related 
to high-frequency gain, but the relationship is not 
obvious. The relationship deduced here is based 
on the damping ratio and peak magnitude of the 
Bode plot. Because the damping ratio C i s related 
to overshoot by 


\\G(jLj)W p (ju;)\U. = 1 (Cl) 


Percent overshoot = 


e-Or/Vl-C 2 


where G{juj) is any transfer function between the 
input and output of interest. From classical control, 
it is known, for example, that if a small steady-state 
error is desired, then the transfer function between 
the error and reference signal should be small in 
the frequency range' of interest. The same is true 
for any set of input/output combinations. Thus, 
if G(ju;) must be small at low frequencies for good 
tracking performance, then, conversely, Wp must be 
large. In addition, from the Bode integral theorem, 
G (, j Us ) must become large at high frequencies, thus 
forcing W p to become small. Based on the criterion 
described in equation (Cl) and -on the relationships 
derived from classical control, the method presented 
in this section can be used to calculate W p . 


Combining specifications of steady-state error, 
percentage of overshoot, ami time constant for the 
variable of interest, the time domain specifications 
can be directly translated into a frequency-dependent 


transfer function in the //^ domain. Common rela- 
tionships from classical control make the time con- 
stant r of the response equivalent to the inverse of the 
crossover frequency for the descriptive transfer func- 
tion. Furthermore, applying the final value theorem 
yields the steady-state tracking error that is analo- 
gous to the inverse of low- frequency system gain. The 
performance weighting function is for a high-pass fil- 
ter. so low- frequency gain can be also interpreted as 
the amount by which disturbances are attenuated. 
Thus, a performance weighting function construction 
proceeds as follows. Assume the transfer function has 
the form 


W p {joj) = K 


, ci juj T 1 


(C2) 


pKj ~' "bju+ 1 
Because low-frequency gain gives a steady-state er- 


ror, let 


0 : K if — K 


1 

C.s.s 


(C3) 


From the Bode diagram at crossover frequency uj c , 


L m \Wp{jv r )\ = 0 implies that \W p {jw r )\ = 1 


Solving for uu r results in 


uv 


1 

r 


1 - K 2 
K 2 ci 2 - b 2 


(C4) 


the peak magnitude of the standard second-order 
transfer function for the desired overshoot can be 
obtained from the Bode plot. As the values of C ap- 
proach critical damping (0.707), the peak magnitude 
approaches 0 dB and this method no longer applies. 
However, for C < 0.707, the inverse of peak magni- 
tude serves as the high-frequency gain /vy, resulting 
in 

Ka 

u: — > oo : Ahf = (C5) 

Hence, equations (C3) (C5) can be solved for the 
three unknowns A', a, and b. Equation (C3) gives the 
value for K. The expression for b from equation (CT>) 
is substituted into equation (C4) and the resulting 
equation is solved for a, yielding 

9 Iruj 2 — IT K 2 

a ~ = WJ}. — 


Thus, 


a = 


K 2 - 1 



Tfr 


(C6) 


For example, to calculate the performance trans- 
fer function for velocity response to velocity com- 
mand, lot 

c.s.s - 5% = 0.05 - K = 20 

r — 40 sec — > u r = 0.025 

10% overshoot — ► £ = 0.6 — > A^f = 0.97 


The resulting transfer function does not guar- 
antee that the variable of interest will follow 7 the 
specified performance. The three response specifi- 
cations are not physically independent of each other 
so even if they are independently designated in the 
performance weighting, the desired response w T ould 
not be attained without knowing the physical capa- 
bilities or inclinations of the variables. In the ex- 
ample problem, the achieved overshoot is a much 
more direct function of the specified time constant 
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than of the specified overshoot in the commanded 
variables. Moreover, for the time constants speci- 
fied in the problem, the commanded variables were 
unlikely to produce large overshoots even when al- 
lowed by the performance weightings (i.e., £ was 
notably relaxed). In fact, the overshoot specifica- 
tion that was incorporated into W p had to be re- 
laxed because it considerably slowed the time con- 
stant. The overshoot incorporated into W p that 
produced the actual specified response was 40 per- 
cent (40% overshoot — ► £ = 0-3 — > = 0.5). Simi- 

lar adjustment had to performed on the steady-state 
error that was incorporated into W p . The desired 
steady-state error was less than 5 percent; the one 
incorporated into W p was 1 percent. 

The time constant for velocity response, remained 
at 40 sec which resulted in 


i,i- i- _ °- 5 (s + 4.330 x nr 2 ) 
p r ~ .s + 2.165 x 10 “ 1 


(C7) 


The altitude response had a propensity for faster rise 
time than was specified. To minimize the elevon 
response that directly affected the angle of attack, 
the altitude time constant was increased to 45 sec, 
which resulted in the transfer function 


0.5(s + 3.849 x 10 - 2 ) 
.s + 1.925 x 10-'* 


(C8) 


The weighting of a was based on the desire to 
attenuate atmospheric turbulence and to indicate the 
import ance of a response in the opt imization process. 


The starting number was 10 and the final constant 
was derived by iteration. However, as mentioned 
above, elevon response has a great deal more effect 
on a response than the actual a weighting. Unity 
weightings on q and 0 indicated those responses are of 
secondary importance. Due to the linear relationship 
linking a, 6 , and 7 (flight path angle) that was 
assumed here, the magnitude of the 9 response was 
similar to that of a. Moreover, q was also linked 
by a linear relation to 9. Consequently, both q and 
9 responses were satisfactory and resulted from the 
constant weightings 


a 


" 20 " 

Q 

= 

1 

_ 9 _ 


1 


(C9) 


The actuator position and rate limits were imposed 

by 


Wp 

1 1 

2. O 

3 ft 

\ 1 

- 

30 

60 

(CIO) 

Wp 

be 
bin j 

- 

20 

10 

(Cll) 


where the position weighting function started as an 
allowable deflection limit for the elevon (± 20 °) and 
for the fuel flow rate (60 percent). To minimize 
elevon deflection, the weighting function had to be in- 
creased. The actuator rate weightings were* similarly 
established. The starting point was the assumed rate 
limit that had to be relaxed for both control variables 
to achieve satisfactory responses. 
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Appendix D 

Numerical Systems 

The longitudinal dynamics are represented by 


x ~ Ax Hb Du 


(Dl) 


where 




a 

Q 

0 

_h_ 


u = 


he 

V 


The system and control matrices are set at a Mach number of 8, an altitude of 85 / 00 ft, and a dynamic pressure 
trajectory of 2000 psf to give 


- 3.65424 x 1CT :} 

-9.6679 x KP 1 

0 

-5.5639 x 10 1 

-1.4321 x 10"' 5 - 

-3.91925 x 10- r > 

-8.1626 x 10“ 2 

1.0000 

-8.4420 x 10~ 5 

9.2560 x 10 -6 

2.0147 x 10- ;J 

3.0354 

-9.5218 x 10~ 2 

1.5500 x 10 5 

-1.0766 x 10 _r> 

2.7263 x 10~ (i 

7.7679 x 10~ {i 

1.0000 

-7.7679 x 10 (i 

-1.0188 x 10 _!) 

. 2.0779 x 10~ 2 

-1.3701 x 10 2 

0 

1.3701 x 10 2 

0 


and 


B = 


-9.6995 x 10“ 2 
3.3486 x 10“ 3 
1.0825 
0 
0 


7.5989 

-2.0942 x 10~ 3 
0 
0 
0 


Aerodynamic forces and moments on the airplane depend on the relative motion of the airplane in the 
atmosphere and not on the inertial velocities. To account for atmospheric disturbances, the forces and moments 
must be related to the relative motion with respect to the atmosphere. This accounting is done by expressing 
velocities used in calculating aerodynamics in terms of the inertial and gust, velocity components. Hence, the 
atmospheric turbulence is introduced conventionally into the state equations as 

x = Ax + Bu + Ev (D2) 


where 


E = 


--3.6524 x 1(T 3 
3.9195 x 1CT 5 
-2.0147 x 10~ 3 
0 
0 


9.6679 x 10“ 1 ‘ 
8.1626 x 10“ 2 
-3.0354 
0 
0 


v — 


e gU st 

a gust 


Here, gust matrix terms 0 and h = 0 because they depend strictly on inertial velocities. 
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The H x controller is a dynamic suboptinml compensator with the number of states // that is (equivalent 
to the number of generalized plant P states; for this problem, n — 12. The state equations represe nting the 
system are 


'(l‘ 




= 






T B ( 


v ■ 

a 

< 1 
0 
h 

Vc 

L h<- 


. . i 

r(.i 

ii 

-c £ 

'C 

zn 

i 


(D3) 


The // controller is based on the // x controller with the generalized plant P scaled by matrices D and D -1 to 
reflect, the structure of the uncertainty and has the same state space form as shown in equation (D3). Thus, 
the new plant P — DPI) 1 is augmented by the states of D and D~ l during each iteration, with the final 
controller containing 18 states. The controller was reduced using the truncation with residue and optimal 
Hankel singular values methods. The resultant 13-state controller produced an essentially identical response 
in a closed- loop system as the original // controller. 

For completeness, the poles of the // x controller are given in the table for the nominal system (referred to as 
the baseline controller), for an // x controller that is designed for 20-percent uncertainty, and for a // controller 
that is based on the latter // -Xi controller. Neither the /f-x- controllers nor the // controller have transmission 
zeros. The closed-loop functions for the controllers are also characterized as having no transmission zeros. 

The // x baseline controller remains robustly stable to 17-percent, uncertainty, so the destabilizing pertur- 
bation in a simulation must satisfy ||U / ^A, S || 0C , = 0.17. The peak value fi[G\ \ (jut)] = 1.174 for 20-percent 
uncertainty implies that the system will remain stable for ||A||rxi = 1/1-17 1 — 0.8b instead of for ||A||-x, = 1. 
This value' also implies that ||l! /^A. s ||x; — ||0.20 • 0.85|| x = 0.17. The // controller not only remains stable but 
also satisfies robust performance for 32. 5-percent, uncertainty (the robust performance criterion, /*[("»( jo;)] 1. 
so the worst ease performance perturbation must, satisfy IIH aA^II^ — 0.325 in a simulation. 
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Nominal System 


Poles 

Frequency 

Damping 

h\(A 

. = 0),n = 12 


2.0407 x 10~' ± j 1.4148 x 10“ ' 

2.0456 x 10 4 

0.998 

1 .3431 x 10“' ± jl.7063 x HP 1 

.21715 

.619 

-2.7473 

-2.7473 

1.0 

-2.5096 ±j A. 9674 

5.5654 

.451 

42.583 

42.583 

-1.0 

— 65.745 ± j'68.746 

95.124 

.691 

-293.56 ± j 275.83 

402.81 

.729 

h\( A 

= 20). n = 12 


-1.9250 x 10 1 

1.9250 x 10“ 4 

1 

-2.1650 x 10 4 

2.1650 x 10~ 4 

1 

-.12428 

.12428 

1 

-.17073 

.17073 

1 

-1.4143 

1.4143 

1 

—6.2865 

6.2865 

1 

-66.508 

66.508 

1 

— 37.734 ± j 60. 294 

71.128 

.531 

-2444.1 

2444.1 

1 

- 1 .0824 x 10° 

1.0824 x 10 fi 

1 

-7.5995 x 10° 

7.5995 x 10° 

1 

tf,i( A = 

= 20 / 2 ), n = 13 


-1.9242 x UP 4 

1.9242 x 10- 4 

1 

-2.1651 x HI' 4 

2.1651 x nr 4 

1 

-8.0876 x 10“ 2 

8.0876 x 10“ 2 

1 

-.17008 

.17008 

1 

-.57388 ± j 0.77693 

.96589 

.594 

-1.4146 

1.4146 

1 

-2.6054 

2.6054 

1 

-8.2357 

8.2357 

1 

-38.874 

38.874 

1 

-20.440 ± j 34. 699 

40.272 

.508 




Appendix E 

Atmospheric Turbulence Model 

The traditional approach to modeling atmo- 
spheric turbulence as a stochastic process has been to 
use the Dryden spectra. For engineering purposes it 
is assumed that the power spectra of atmospheric tur- 
bulence can be approximated by the Dryden spectra. 
( See ref. 1(5.) The longitudinal (body-axis) Dryden 
power spectrum is defined by 




Lu_ 2rrJ 1 

\ 7r 1 + (L u y) 2 


and the lateral and vertical is defined by 


(El) 


Lr.w Vr.w- 1 +3(L r . jr ^)- 

>' ~fi+UW)*|* {h2> 

Here, V is the' vehicle 1 velocity (unit length per sec- 
ond) and uj is the spatial frequency and is in radians 


per second. I he variance of the turbulence is a 2 
and the spatial scale* length is L. Increasing a scales 
the ])o\ver spectral density (PSD) or if increases the 1 
poweT at all freepuncies without changing the 1 rela- 
tive distribution. Changing L redistributes the poweT 
over the frequency range 1 . Increasing the 1 scale 1 length 
of the 1 turbulence* incTeases the ]>ower at lower fre- 
quencies and ile*croases the power at higher freepien- 
cies such that the 1 integral of the 1 PSD remains con- 


stant; the 1 integral of the PSD s will alwavs expiate 1 to 

) r . 1 

rc tor /('re)- mean input. 


The 1 exponents of the Dryden spectra are 1 int.eg(*rs, 
so filteTs can be 1 developed through which unit vari- 
ance 1 Gaussian white noise may be passexl. The 1 re*- 
sulting filtere'd white noise secpiencc* will have 1 ne'arlv 
the same 1 statistical pniperties as re*corded turbulence 
time historic's. 


To determine the 1 filter expiation. the following 
relation is used: 


= |//(»| 2 *in (E3) 

wIutc ‘Font is the elesireel Dryden spe'ctrmn, H (ju) 
is the 1 unknown filter or transfer function, and <I>j n 
is the 1 poweT spexlrum of the input white noise 1 se 1 - 
epiemce*. In eontinuous time 1 , white 1 noise* contains all 
frexpiencie's so the integral of the PSD is not defineel. 
The magnitude of the PSD across all frequencies is 
prescribed t o be unity. Therefore, in continuous time 
applications, the 1 filters are simply the spectral roots 
of the Dryden spex-tnnn expiations or 


Hu (j&) — 


1 

J — f j- 


(E4) 


for the longitudinal turbulence 1 and 


:iV 


Hr.tr(j^) — \j 

y I'rjc'K 


(S“~i 




(E5) 


for the 1 lateral and vertical turbulence. 

In the discrete. 1 time domain, as in a digital simu- 
lation. -pure, white noise cannot be geneTatexl to pass 
through filters; the simulation can generate frexpien- 
e‘ie\s only up to the Nyepiist frexpnnev. The* result is a 
banel-limited white noise sequence* with a PSD that 
is constant but not unity. The 1 intensity is now deter- 
miiK'el from the requirement that the 1 integral must 
be 1 the 1 variance that is usually pn*se*ribe*ef to be 1 unity. 
Therefore. the intensity e>f the 1 input noise* sexpience 1 
is the inverse of the Nyepiist frequency or (Ai)/n. 

The 1 values of a anel L are tabulatex! in numerous 
references as functions of altitude anel longitudinal, 
ve'rtical. and lateral components. For (he 1 r( , f(*r<*nce 
altitude usexl in this stuely (85 700 ft), the <r*s anel 
//s for lemgit uelinal and vertical turbulence (red. 17) 
(lateral was not modelexl) we're 1 


er u — 10.8 ft/sev rr (r = (i.88 ft /se*c 

L (l = 05 574 ft L tr - 20 220 ft 


With the'se 1 filter parame'teTS implementexl in the 1 
Laplace* domain for a discrete time 1 simulation, the* 
longitudinal turbulence* filte*r bexome's 




/X /X'v 

V Af y /.|/7T 



(K«) 


and F (r dccoiiK's 


/V(.s) = 


/X / :U \ :i 4 • s ) 

\ At\j L„, 7 T (y+.s)^ 


(K7) 


The use of filt.erexl white 1 noise 1 to represent t urbulonce* 
to te'st or evaluate 1 a design make's verification of ce>r- 
re'c-t ly geneTatexl sirnulatexl turbulence 1 ne'cessary. For 
a time 1 simulat ion, theemtput turbulence time* history 
should be aimlvzexl and the* standard deviation of the 
signal should equal (to a dost* approximation for a 
sufficient, number of time* samplers) t he value* usexl in 
the filter. Furthermore 1 , the 1 integral e)f any PSD of 
the filter output, theoretical or experimental, should 
equate to a 2 . 
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modt l. an integrated flight propulsion dynamic model of a conical accelerator vehicle, was linearized as the vehicle accelerated through a Mach 
number of S. Vehicle acceleration through the selected flight condition gives list' to parametric variation that was modeled as a structured 
i meet t aint.v. I hr p analysis apprinuh was used in the frequency domain to conduct controller analysis and w r as confirmed by time history plots. 

I In' results from (his example demonstrate the inherent advantages of I he p framework for t his class of problems. 
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